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SUMMARY 


The empennage structure on the Learjet 55 aircraft was quite similar to the 
empennage structure on earlier Learjet models. However, due to an important struc- 
tural change in the vertical fin along with the new loads environment on the 50 
series aircraft, a structural test was required on the vertical fin, but the hori- 
zontal tail was substantiated by a comparative analysis with previous tests. 

NASTRAN analysis was us»d to investigate empennage deflections, stress levels, and 
control surface hinge forces. The hinge force calculations were made with the 
control surfaces in the deflected as well as undeflected configurations. A skin 
panel buckling analysis was also performed, and the non-linear effects of buckling 
were simulated in the NASTRAN model to more accurately define internal loads and 
stress levels. Comparisons were then made between the Model 55 and the Model 35/36 
stresses and internal forces to determine which components were qualified by previ- 
ous tests. Some of the methods and techniques used in this analysis have been 
described in the following paragraphs. 


INTRODUCTION 


Features in the 55 empennage 'section that were identical to previous Learjet 
models included the horizontal tail, elevators, and the upper portion of the verti- 
cal fin and rudder (see figure 1). The most significant change in the empennage 
occurred in the lower portion of the vertical fin and rudder which were expanded 
using the same taper ratio and airfoil contour as the upper portion so that the 
total height of the vertical fin and rudder was increased by approximately 15 inches 
Since the vertical tail increased in size and the vertical fin spar depths and 
spacings changed at the intersection with the tailcone, an entirely new tailcone 
section was designed to support the empennage section. Control systems for the e^e 
vators and rudder were nearly the same as the previous controls in this area with 
the exception of the 15 inch height increase in the vertical tail and the geometry 
changes in the tailcone. 

Construction in the vertical tail consisted of five spars and eight ribs made 
of bent up aluminum sheet metal covered with an aluminum skin. The rudder was 
assembled with one spar, eleven ribs, and a thin skin cover, all of which were alum 
inum sheet metal fabrication. Support for the rudder was accomplished through 
three hinge points. Two of these hinges, the top and center attachments, connected 
the rudder leading edge to the rear spar of the vertical stabilizer while the third 
hinge point, or bottom support, was attached to the aft side of frame 48 through a 
torque tube and bracket support assembly (see figure 2). Near the lower end of 
this torque tube was mounted a bellcrank assembly which provided the means for 
actuating and controlling the rudder surface. 
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structural arrangement of the horizontal tail included a forward and rear spar 
with three intermediate stringers on both the upper and lower surfaces, seven ribs 
per side, and upper and lower skin panels. The elevators were one spar airfoil 
sections with nine ribs per side which were covered with an aluminum skin. Each 
elevator was attached to the horizontal stabilizer by three hinge points along the 
elevator leading edge. The uip and center span iiinges both connected directly to 
the rear spar of the horizontal stabilizer, but the inboard hinge attachment was 
very similar to the lower hinge on the rudder. 'I'he elevator inboard hinge consisted 
of a torque tube assembly which was supported at the horizontal stab'lizer rear spar 
and aircraft center line by a bracket assembly. This bracket installation served 
as a support for bo'^h the right hand and left hand elevators at the inboard end 
(see figure 3). Actuation and control of the elevators was achieved through indi- 
vidual bellcranks on the inboard end of each torque tube. These bellcrank asserri- 
blies were then connected by a push rod linkage system which extended through the 
vertical stabilizer to a control system sector in the fuselage tailcone. 


SYMBOLS 

A-j X, Y, and Z components of a point at the origin for a C0RD2R card 

a rotation about the x axis for element internal forces 

Bi X., Y, and Z components of a point on the Z axis for a C0RD2R card 

e rotation about the Y axis for element internal forces 

Ci X, Y, and Z components of a point in the X-Z plane for a C0RD2R card 
internal force component 

Y rotation about z axis for element internal forces 

Hq coordinate values of a deflected control surface transformed to the 
basic system 

coordinate values of a deflected control surface in a local system 
coordinate values corrected for hinge line sweep 
He coordinate values corrected for control surface rotation 
A hinge line sweep 

Ta transformation matrix for hinge line sweep 

T0 transformation matrix for control surface rotation 

e control surface rotation 

X first coordinate of a point in a local rectangular system 

y second coordinate of a point in a local rectangular system 
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third coordinate of a point in a local rectangular system 
X first coordinate of a point in the basic rectangular system 

Y second coordinate of a point in the basic rectangular system 

Z third coordinate of a point in the basic rectangular system 


BACKGROUND 


Since the Model 55 empennage structure was very similar to the Model 35/36 
empennage structure and other earlier empennage structures, Learjet was faced with 
the situation on the 50 series development program of determining which structural 
items were qualified by previous certification tests and which structural items 
needed to be requalified, .’he vertical tail and tailcone support structure was 
obviously one area which had to be qualified by structural test, since the height 
of this component had been increased on the 50 series aircraft. A portion of the 
NASTRAN analysis on that structure has been detailed in NASA CP-2151 (see ref. 1) 
and consequently has not been addressed in this :ext. However, a description of 
the NASTRAN analysis performed on the rudder and rudder hinges was not covered in 
that document and has been incorporated in the body of this paper. 

A review of the loads for the 50 series aircraft helped to further clarify 
some of these questions, but for some structural components the answer was not 
obvious and more study was necessary. One area where this was particularly true 
was the horizontal tail installation. This situation was mainly due to a revised 
loads definition for the horizontal tail where the combination of new loads was not 
directly comparable to the older load conditions. Also, a faster method was needed 
for calculating elevator hinge forces for the various load conditions and deflected 
elevator configurations. Previously, hinge forces had been calculated using methods 
which were very tedious, time consuming, and not very adaptable to changes in 
elevator position. 

To help clarify these questions a NASTRAN analysis was performed on the empen- 
nage structure for the various load conditions and control surface positions so 
that a comparison could be made between the deflections, stresses, and internal 
loads on the 55 empennage with the 35/36 empennage. Previous structural tests had 
demonstrated that the horizontal stabilizer and elevators on the 35/36 aircraft 
had the capability to sustain 150% of ultimate load (not limit load) for two of 
the more critical conditions on that airplane. The ib0% of ultimate load tests 
were performed to more accurately define the strength of the horizontal stabilizer 
and elevators as well as establish the growth capability of this structure. Con- 
sequently, a portion of the NASTRAN analysis was spent simulating the 35/36 load 
conditions in order to establish a baseline for the 50 series analysis. Then the 
50 series loads were applied to the NASTRAN model and a comparison was made between 
these data and the 35/36 data. Using this approach Learjet was able to more 
clearly define which elements were qualified by previous certification tests and 
which members, if any, need to be qualified in the 50 series substantiation pro- 
gram. This data base would also help to establish what components needed to be 
qualified in future growth versions and/or modifications to the existing config- 
uration. 
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MODELING CRITERIA 


The NASTRAN model used to represent the empennage section included the ver- 
tical stabilizer, rudder, horizontal stabilizer, elevators, and tailcone section 
(see figure 1). The tailcone was utilized in the model to more accurately des- 
cribe the flexibility and stiffness of the vertical tail, since the spars from the 
vertical tail extended into and were attached to the tailcone. A significant por- 
tion of the control system connecting the rudder and elevators was also modeled 
to more accurately define the effect of control system flexibility on the deflec- 
tion of the control surfaces and hinge forces. 

Since this structure was geometrically quite complex, multiple local roor- 
dinate systems were used to simplify the modeling and to make the analysis easier 
to perform for the different control system deflected configurations. Each frame 
in the tailcone section was modeled in a local cylindrical coordinate system while 
the vertical fin and horizontal stabilizer were modeled in a local rectangular 
coordinate system. The rudder and elevators were all modeled in local rectangular 
systems with the Y axis of each coordinate system defined along the hinge line. 
This type of coordinate system definition permitted the rotation of each control 
surface about the hinge line by simply transforming the C0RD2R card (see ref. 2) 
rather than changing all of the grid point coordinates in the control surface to 
reflect the deflected position. Rotation of the control surface about the hinge 
line thus permitted the calculation of hinge forces for various control surface 
movements. This technique also facilitated a study of vertical fin and horizon- 
tal stabilizer internal forces and stresses for various control surface positions. 

The model for this structure was approximately 5700 degrees of freedom in 
size with 950 nodes and 2800 elements. Local coordinate systems totaled seven- 
teen which consisted of both rectangular and cylindrical reference systems. Since 
the results of previous static tests had indicated that skin panel buckling was 
present, a buckling analysis was implemented to identify buckled skin panels and 
to simulate this non-linear phenomtcon in the NASTRAN analysis (see ref. 3). This 
buckling simulation was achieved by replacing the original skin thickness of the 
buckled panel with an effective skin thickness which represented the stiffness of 
the buckled configuration. The buckling analysis and effective skin thickness 
calculation were performed by a computer program which was set up to run as a post 
processor to the NASTRAN program. This process generally had to be repeated 
several times in order to obtain a convergent solution. A more detailed descrip- 
tion of the theory behind this analysis has been presented in NASA TMX 3428 and 
further discussion will not be covered in this text. 


TAILCONE 


Incorporation of the tailcone section with the empennage model was done to 
more accurately define the stiffness relationship between these two components. 
Representation of the tailcone section included the st.ucture from frame 39 to 
frame 48 (see figure 1). The modeling of this structure was accomplished by 
defining the grid points at each frame and stringer intersection in a local cylin- 
drical coordinate system in the plane of the frame using C0RD2C cards (see ref. 2). 
NASTRAN members used to model this structure were BAR elements for the frames, 
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doublers. A buggaye uuui oi, 1.1, n. 11 , iu siue ut ihe tailcone did not have sig- 
nificant structural carry through capability and was omitted to conserve degrees 
of freedom. The keel beam structure that extended from frame 39 to 41 in this 
section was modeled using CONROD elements for the caps, SHEAR eVments for the side 
webs, and QDMEM2 elements •^or the horizontal webs. 


VERTICmu tail 


Two separate local rectangular cocrdi.iate systems were used to facilitate the 
definition of grid point locations in the vertical fin and rudder (see figrre 2). 

The local rectangular coordinate system for the vertical fin was defined using a 
C0RD2R card with the origin located on the center line of the airplane at the inter- 
section of the vertical tail leading edge and the tailcone upper surface. The X 
axis was defined as positive aft, the Y axis was defined as positive up, and the Z 
axis was defined as positive left hand outboard. The local rectangular coordinate 
system for the rudder was defined using a C0RD2R caro with the origin on the rudder 
hinge line at the lower end of the rudder torque tube. Orientation of the '' axis 
was positive aft and slightly down (for zero degree rudder deflection), the Y axis 
was positive up along the rudder hinge line, and the Z axis was positive left hand 
outboard. 

Grid points were located at the intersection of the spar caps and rib caps 
on both the right hand and lef-^ hand outer contours of the vertical tail. The 
spar caps in the region of the fuselage frames were modeled using BAR elements 
while the remainder of the spar caps and rib caps were simulated by CONROD elements. 
Skin panels were represented by QDMEM2 members, although there were some triangu- 
lar, TRMEM, membrane elements used to describe intricate corners. SHEAR elements 
were used to represent .;he spar and rib webs. 

Attachment of the rudder to the vertical fin was accomplished through two 
hinges from the aft spar of the vertical fin to the leading edge of the rudder. 

The hinge model consisted of a framework of eight CONRODs that connected the hinge 
point to a rudder rib on the aft side and to a vertical fin rib on the forward 
side. The first group of four CONRODs, or the aft hinge, connected the hinge point 
to four points on the adjacent rudder rib, while the second group of four l,ONRODs, 
or the forward hinge, connected the hinge point to four grid pc ints cn the adjacent 
vertical fin rib (see figure 2). This method of modeling facilitated an easier 
computation of hinge loads through a resolution of forces in the CONRODs into 
forces in a predefined coordinate system by using a series of transformation equa- 
tions. The torque load in the rudder was restrained by a bellcrank mounted on a 
torque tube whose axis coincided with the hinge line of the rudder. This torque 
tube was attached to the bottom of the rudder on one end and to a bracket support 
on the aft side of frame 48 on the other end. BAR elements were used to model the 
torque tubes and the bellcrank assembly which actuated and controlled the rudder 
rotation at the bottom of the torque tube. Modeling of the bracket support on 
frame 48 was accomplished using CONRODs for the stiffeners and QDMEM2 membranes 
for the webs. 




230 


\ny.<: 


All. 


The horizontal tail consisted of the horizontal stabilizer and right hand and 
left hand elevators. Attachment of the horizontal stabilizer to the vertical fin 
was accomplished through a hinge point on the top of the vertical fin between spars 
three and four and a screwjack attached to the horizontal stabilizer apex and spar 
two of the vertical fin. Rotation of the horizontal tail could then be achieved 
through actuation of this screwjack. The elevators were attached to the rear spar 
of the horizontal stabilizer through three hinges on each side. 

Three separate local rectangular coordinate systems were created using the 
C0RD2R facility in NASTRAN to model the horizontal tail structure. These local 
coordinate systems were used to describe the horizontal stabilizer, right hand ele- 
vator, and left hand elevator. Layout of the horizontal stabilizer local rectan- 
gular coordinate systems was established with the origin located at the horizontal 
stabilizer apex with the X axis positive aft, the Y axis positive left hard out- 
board, and the Z axis positive down. Orientation of the left hand elevatoi local 
rectangular coordinate system was defined with the origin at the intersection of 
the hinge line and aircraft center plane with the X axis perpendicular to the hinge 
line and positive a^t, the Y axis coincident with the hinge line and positive left 
hand outboard, and the Z axis perpendicular to the X-Y plane and positive down. The 
right hand elevator local rectangular coordinate system was similar to the left hand 
elevator with some differences in positive axes directions. The origin for the 
right hand elevator coordinate system was the same as the left hand elevator; how- 
ever, the X axis was perpendicular to the right hand hinge line and positive for- 
ward, the Y axis coincided with the right hand hinge line and was positive right 
hand outboard, and the Z axis was perpendicular to the X-Y plane and was positive 
up. One important advantage of the separate coordinate systems for the elevators 
was the ease and improved flexibility of analyzing the elevators and horizontal 
stabilizer with the elevators in various deflected positions. 

Grid points in the horizontal tail were specified in relation to the appro- 
priate coordinate system depending on whether the point v.:s in the horizontal 
stabilizer, “"ight hand elevator, or left hand elevator, .ocation of ihe grid 
points was defined by the intersection of the spars and stringers with the ribs 
and the exterior contour. As in the vertical fin, the spar caps, stringers and 
rib caps were represented by CONRODs, the spar and rio webs were SHEAR elemeiits 
and the skin panels were modeled using QDMEM2 and TRMEM members in the horizontal 
stabili/ier while QUAD2 and TRIA2 elements were used in the elevators. BAR ele- 
ments were used to simulate the elevator trailing edges. 

Attachment of the elevators to the horizontal stabil' ar rear spar was accom- 
plished by three hinges on each elevator. These hinges were designated by the 
spanwise location and were referred to as the inboard hinge, center hinge, and tip 
hinge (see figure 3). The center and tip hinges were modeled as a framev/ork of 
eight CONRODs that connected the hinge point to an elevator rib and a horizontal 
stabilizer rib. The first group of four CONRODs, or aft hinge, connected the hinge 
point to four grid points on an elevator rib, while the second group of four 
CONROns. nr forward hinge, connected the hinge point to fGur grid points on a 
horizontal stabilizer rib (see figure 3). The torque load on the elevator was 
reacted by a belicrank mounted on a torque tube whose axis coincided with the ele- 
vator hinge line. There were two of these bel Icrcnk/torque tube assemblies, one 
for each side. The torque tube attached to the elevator on the tube's outboard 
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end and to a bracket support mounted on the horizontal stabilizer's rear spar on the 
airplane center line on the tube's inboard end. Both torque tubes and bellcranks 
were modeled using BAR elements, while the center line bracket support assembly 
was modeled using CONRODs. 


CONTROL SYSTEM 


A portion of the control system connecting the rudder and elevator has also 
been included in the modeling of the empennage section to more accurately define 
the effect of control system stiffness on the deflection of these control surfaces 
and in particular the associated hinge forces. Rudder motioii on the Learjet M55 
aircraft was accomplished by rudder pedal movement which was transferred by a 
closed loop cable assembly which drove a bellcrank mounted on the rudder torque 
tube (see figure 4). The aft portion of this cable system was included in the 
finite element moael and was represented by CONROD elements with the side of the 
closed loop cable system that was not in tension removed. The CONROD elements 
representing the cables were extended into the tailcone section where the forward 
end of this member was attached to a pulley which was restrained against rotation. 

Elevator deflection on this aircraft was achieved by control column movement 
which was transferred through a closed loop cable assembly and a push-pull rod 
linkage system. The upper portion of the push rod system consisted of two round 
tubes which were attached to the bellcranks on the left and right elevators respec- 
tively. The other end of these two push roas were joined to a down spring linkage 
mounted on the aft side of spar three in the vertical fin. A single round tube was 
used in the lower portion of the push rod system. This shaft was connected to the 
down spring linkage on the top end and to a sector on the bottom end (see figure 5). 
The portion of the elevator control system that was included in the empennage 
model consisted of the elevator bellcranks, upper push rods, down spring linkage, 
and lower push rod which was constrained at the lower end for the NASTRAN analysis. 
The grid points that defined the top and bottom points of the lower push rod were 
defined in separate local cylindrical coordinate systems. This was done so that 
these grid points could be moved with the elevator motion by using a simple set of 
transformation equations. The push-pull rods were modeled with CONRODs, whereas 
the bellcranks and down spring linkage were represented by an assemblage of BAR 
elements. 

The horizontal stabilizer rotation about the hinge point between spar three 
and four on top of the vertical fin was achieved by a screwjack actuator attached 
to the horizontal stabilizer apex. A BAR element was used to represent the screw- 
jack. One end of this BAR element was attached to a fitting on the forward side 
of spar two in the vertical fin, and the other end was connected to a fitting 
mounted on the horizontal stabilizer center line. The fittings that attached the 
screwjack actuator at both ends were also modeled using BAR elements. 
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LOADS AND CONSTRAINTS 


Many different load conditions were analyzed using the NASTRAN empennage model. 
A correlation analysis was performed on the horizontal tail using a positive gust 
with top roll moment load condition from the 35/36 fatigue test spectrum. This load 
case was selected basically because earlier tests did not have sufficient instru- 
mentation to perform a correlation analysis, and this particular case was one of 
the higher load conditions in the 35/36 fatigue test spectrum. Since the corre- 
lation analysis for the vertical tail was described in a previous NASA paper (See 
ref. 1), the results of that work were not presented here. The analysis of the 
Learjet 55 vertical tail described in this paper was directed mainly at rudder 
hinge force calculations. The load case applied to the vertical tail for these 
rudder hinge force calculations was a maximum side bending with rudder kick load 
condi tion. 

NASTRAN runs were made for all of the 35/36 horizontal tail static test loads. 
These conditions were run in order to establish a baseline for the Model 55 analy- 
sis so that a determination could be made as to whether any of the 55 horizontal 
tail internal loads and stresses exceeded those on the 35/36 horizontal tail. Since 
some of the 35/36 load cases were tested to 150% ultimate load, there was good 
reason to believe that the 55 horizontal tail load conditions were covered by the 
previous 35/36 static tests. The 35/36 loads were applied to the horizontal sta- 
bilizer and elevators with the elevators in the neutral, or undeflected position, 
the 15° up position and the 15° down position. The loads definition for the Model 
55 horizontal tail was somewhat different from the 35/36 horizontal tail. There 
were many more load cases for the Model 55 airplane than for the 35/36 airplane, 
and a direct comparison was not easily achievable for many of these conditions. 
Consequently, internal loads and stresses were determined for the 55 horizontal tail 
load conditions with the elevators in the undeflected and deflected positions using 
NASTRAN analysis in order to establish a common basis for comparison with the 35/36 
internal loads and stresses. 

Loads were applied to the vertical fin, rudder, horizontal stabilizer and 
elevators using FORCE card images (See ref. 2). These forces were applied at the 
grid points to simulate the shear, bending and torque envelopes for these air- 
foil surfaces. A computer program was written which accepted concentrated forces 
from the test loads data as input and redistributed these forces proportionately 
to the four nearest grid points according to the distance from those nodes and 
generated FORCE card images as output. 

Constijints for the empennage model were applied at several locations. The 
majority of constraints were applied along the perimeter of frame 39 (See figure 2) 
in all three translational degrees of freedom at the frame and stringer inter- 
sections. Other constraints were applied to the rudder control system and elevator 
^''ntrol system. The rudder control system was constrained by attaching the forward 
end of the tailcone rudder cable to a frame in the tailcone section. The elevator 
control system was restrained by using SPC, single point constraints, (See ref. 2) 
on the three translational degrees of freedom at the bottom end of the lower push 
rod (See figure 2). 
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Grid coordinates on the control surfaces were specified with respect to the 
appropriate rudder or elevator local rectangular coordinate system. These local 
coordinate systems were initially defined for the control surface in the undeflected 
position. However, the control surface could be deflected or rotated without 
changing the grid point coordinates merely by rotating the local coordinate system. 
This was accomplished by modifying the C0RD2R card image in the bulk data deck by 
transforming the original axes reference to a new axes reference. 

The C0RD2R card image in the bulk data deck for each control surface defined 
a rectangular coordinate system by reference to three points. The first point "A" 

(Al, A2, A3) defined the origin of the local coordinate system (See figure 6). The 
second point "B" (Bl, B2, B3) defined the direction of the "Z" axis, and the third 
point "C" (Cl, C2, C3) defined a point in the X-Z plane. The basic coordinate system 
was used as the reference system for these points (See ref. 2). The origin for the 
basic system was defined at fuselage station (F.S.) 0.00, butt line (B.L.) 0.00, 
and waterline (W.L.) 0.00 with the "X" axis positive aft, "Y" axis positive left 
hand outboard, and the "1" axis positive down. Components of A, B, and C were defined 
on the C0RD2R card image by means of a simple series of transformations described by 
the following matrix equations. 

[%] - M’' K] 

The equation shown above converts coordinates in a swept deflected control 
surface to coordinates in a swept undeflected control surface. Equation (2) shown 
below transformed the coordinates in equation (1) to coordinates parallel to a 
fuselage reference system corrected for the hinge line sweep angle X. 

N - [«e] 

Substituting equation (1) into equation (2) and then referencing these values 
to the basic system with an matrix resulted in equation (3). 

[Hg] = [Tl] ^ [^b] ^ [»l] ^ [Hgo] 

Equation (3) was used to calculate the components of points 
C0RD2R card which defined the rotated control surface. 


(3) 

"B" and "C" on the 


HINGE LOAD CALCULATION 


As indicated previously, the hinges were modeled using CONRODS and/or Bar elements. 
Since the hinge point had to be in equilibrium under the action of the internal element 
forces, the resultant of the forces from the elements on the aft side of the hinge line 
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must be equal and opposite to the resultant from the elements on the forward side 
of the hinge line. This force on the forward hinge member was defined as the hinge 
force for each of the control surface support points. 

Internal forces In the CONRODs and BARs were extracted from the NASTRAN output 
for a given load case. The forces In the hinge elements were then rotated into a 
common coordinate system through a series of transformations (See figure 7). The 
three force components in the element coordinate system were then defined as FxEi Fy£, 
and FzE- Orientation of the element coordinate system was defined by the two end 
points of the member and a ^ vector, if the element was a BAR (See ref. 2). Ele- 
ment grid point components x/\, y/\, Zp^ and xg, yg, zg (See figure 7) established the 
basis for transforming the element forces to a common coordinate system. This was 
accomplished by rotating, or transforming, these forces through three angles a, 3. 
and Y* Upon calculation of these three rotations, the transformation matrices were 
formed as shown in figure 7. The element forces F^r* Fy£, and Fze were then resolved 
in the cjmmon X, Y, Z directions, and then the hinge force components were determined 
hy <ui.ining the resolved individual element forces. 

A small comp.iter program was written to perform the transformation from the 
element coordinate system to the common reference coordinate system. This program 
also accounted for the effect of control surface rotation. Hinge loads were com- 
puted on both sides of the hinge line as a check on the transformation calculations. 
The hinge loads at all the hinge points plus the loads at the intersection of the 
torque tubes and bellcranks were also summed to check the equilibrium of these hinge 
loads with the control surface applied loads. 


BUCKLING ANALYSIS 


A comparison of skin cover compressive stresses with the corresponding panel 
buckling allowables indicated the presence of skin panel buckling for several load 
cases. Consequently, an automated buckling analysis was implemented to identify 
buckled skin panels and to simulate this non-linear phenomenon in the NASTRAN analy- 
sis. The buckling simulation was achieved by replacing the original skin thickness 
of the buckled panel with an effective skin thickness which represented the stiff- 
ness of the buckled configuration. A detail description of the theory for this 
analysis has been described in reference 3. The buckling analysis and effective 
skin thickness calculation were coded into a computer program which was set up 
to run as a post processor to the NASTRAN program. This NASTRAN analysis with the 
buckled effective skin thickness generally had to be repeated several times in order 
to obtain a convergent solution. Since some of the loading cases on the horizontal 
tail were combined conditions, the buckling analysis was performed on both the upper 
and lower surfaces simultaneously. 


ANALYTICAL RESULTS 


After extensive analysis and study, the applied loads on the Model 55 hori- 
zontal tail were determined to be less than the maximum tested loads on the 35/36 
horizontal tail which was identical to the 55 in structural detail. Consequently, 
no static test was performed on this structure during the 55 testing program. 
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However, in order to establish a correlation between the horizontal tail portion of 
the NASTRAN empennage model and the structural test data, a positive gust condition 
with top roll moment load case was used from the Model 35/36 fatigue test program 
as was mentioned previously. The results of this test along with the NASTRAN data 
have been exhibited in figures 8 through 11. Figures 8 and 9 show the correlation 
on the horizontal stabilizer forward spar upper cap and lower cap respectively 
while figures 10 and 11 show the correlation on the rear spar upper cap and lower 
cap respectively. Most of the experimental points agree with the NASTRAN data very 
closely with the exception of the most inboard point on the rear spar lower cap. 

This gage was very near an access door in the skin panel and appeared to be affected 
by the stress gradient due to this cutout. As mentioned earlier, no correlation 
was presented on the vertical tail, since this was performed in an earlier paper {See 
ref. 1). 


CONCLUDING REMARKS 


A method has been demonstrated for calculating control surface hinge forces 
using a Learjet Model 55 empennage NASTRAN analysis. Hinge element internal forces 
were extracted from the NASTRAN analysis and converted to a common coordinate system 
using a set of transformation equations to define the individual hinge forces. 

Since each control surface was defined in a local rectangular coordinate system, 
rotation of the control surface could be easily accomplished by transforming the 
components on the C0RD2R card. This control surface rotation capability permitted 
Learjet to investigate the impact of control surface movement on the hinge attach- 
ments in more detail than was previously possible. A buckling analysis was also 
performed to determine the non-linear effect of skin panel buckling on stringer and 
spar cap stresses. These techniques proved to be a valuable asset to Learjet in 
studying the structural operating characteristics of the Model 55 empennage struc- 
ture and associated control surfaces and for substantiating the horizontal tail 
structure by comparison with previous tests. 
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FIGURE 1 - NASTRAN 55 EMPENNAGE MODEL 
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FIGURE 3 - HORIZONTAL TAIL STRUCTURE i 
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Input Data Card CgRD2R 


BULK DATA DECK 
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Rectangular Coordinate System Definition 


Description ; Defines a rectangular coordinate system by reference to the coordinates cf three 
points. TRe first point defines the origin. The second point defines the direction of the 2 -axls. 
The third point defines a vector which, with the z-axls, defines the x-z plane. The reference 
coordinate must be Independently defined. 




F ield Contents 

CIO Coordinate system Identification number (Integer > 0) 

RID Reference to a coordinate system which Is defined Independently of new coordin- 

ate system (Integers 0 or blank) 

A1.A2.A3 

B1,B2,B3 Coordinates of three points In coordinate system defined In field 3 (Real) 

C1.C2.C3 


Remarks ; 1. Continuation card must be present. 

2. The three points (A1, A2, A3), (B1, B2, B3). (Cf , C2, C3) must be unique and non- 
col linear. Noncol linearity Is checked by the geometry processor. 

3. Coordinate system Identification numbers on all C0RD1R, CgRDIC. C0RD1S, CCRD2R. 
C0KD2C, and CgRD2S cards must aM be unique. 


4. An RID of zero references the basic coordinate system. 

5. The location of a grid point (P In the sketch) In this coordinate systt-- cui-n 
(X, V. Z). 

6. The displacement coordinate directions at P are shown py (u^, u^, 

FIGURE 6 - C0RD2R COORDINATE SYSTEM CARD IMAGE 
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FIGURE 7 - HINGE FORCE TRANSFORf-IATION 
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FIGURE 9 - HORIZONTAL STABILIZER FORWARD SPAR LOWER CAP STRESSES 
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HORIZONTAL STABILIZER POSITIVE GUST WITH TOP ROLL MOMENT LOAD 
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AXIAL 

COMPRESSION 

STRESS 
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HORIZONTAL STABILIZER POSITIVE GUST WITH TOP ROLL MOMENT LOAD 

Horizontal Stabilizer Fatigue Test Stress Values, 
Horizontal Stabilizer NASTRAN Stress Values. 



FIGURE 11 - HORIZONTAL STABILIZER REAR SPAR LOWER CAP STRESSES 
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